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SUMMARY 

Algebraic grid generation methods based on transfinite interpolation 
called the two-boundary and four-boundary methods are applied for generating 
grids with highly complex boundaries and yields grid point distributions that 
allow for accurate application to regions of sharp gradients in the physical 
domain or time-dependent problems with small length scale phenomena. Algebraic 
grids are derived using the two-boundary and four-boundary methods for applica- 
tions in both two- and three-dimensional domains. Grids are developed for dis- 
tinctly different geometrical problems and the two-boundary and four-boundary 
methods are demonstrated to be applicable to a wide class of geometries. 


INTRODUCTION 

Many different techniques have been developed for generating computational 
grids required in the finite difference or finite element solutions of partial 
differential equations on arbitrary regions. An emerging problem, however, is 
the generation of grid systems on which solutions can be obtained for complex 
boundary geometries. The importance of the choice of the grid is well known. 

A poorly chosen grid may cause results to be erroneous or may fail to reveal 
critical aspects of the true solution. Any partial differential equation 
(e.g., the Navier-Stokes equations) expressed in a Cartesian coordinate system 
can be transformed to a uniform grid in computational coordinate system, called 
the boundary-fitted coordinate system. The objective of the grid generation 
method is to provide the Jacobian matrix describing the transformation. Alge- 
braic grid generation methods are highly advantageous for direct computation of 
the physical domain as a function of uniform computational grid. These methods 
generate boundary-fitted coordinate systems by algebraically defining distinct 
boundaries of the physical domain and interpolating between these boundaries 
(1) to (6). The algebraic boundary curves can be analytical function or numer- 
ical interpolating functions. These algebraic grid generation methods can gen- 
erate boundary-fitted coordinate systems without the need to solve the partial 
differentia! equations. Precise control of the distribution of grid points in 
the physical spatial domain can be achieved by stretching functions (7) to (9). 
The Jacobian matrix of the transformation can be obtained rapidly and effi- 
ciently by direct analytical differentiation or by numerical differentiation. 
This is in contrast to the di fferenti al -equation grid generation method where a 
partial differential /system is numerically solved by iterative techniques for 
the mapping between the coordinate system of the physical spatial domain and 
the boundary-fitted coordinate system. The transformation matrix (the Jacobian 
matrix) must be obtained by numerical differentiation, and a grid change 
requires a new solution of the partial differential equation system (10) to 
(12). This approach can be very expensive for problems with three-dimensional 


and time-dependent boundary-fitted coordinate systems. Once the boundary- 
fitted coordinate system has been generated by an algebraic or a differential- 
equation grid generation method, the set of governing equations, such as the 
Navier-Stokes equations, are transformed analytically via a generalized coordi- 
nate transformation. This results in a set of governing equations in which 
all derivatives are with respect to the boundary-fitted coordinates. For the 
finite-difference method of solution, the derivatives in the transformed 
governing equations are replaced by finite-difference formulas (13) to (15) to 
form the desired system of f ini te-difference equations. 

Algebraic grid generation method based on transfinite interpolation called 
the two-boundary and four-boundary methods are used to generate grids for two 
and three-dimensional spatial domains with complex arbitrary geometries. A 
detailed description of the methods is presented in reference 1. These methods 
are highly versatile and have a wide variety of applications in both two and 
three dimensions. In the two-boundary method, two separate nonintersecting 
boundaries of the physical domain are defined by means of algebraic functions 
(or numerical interpolation functions). These functions have as independent 
variables, coordinates which are normalized to unity. The two boundaries that 
must be mapped correctly from the physical spatial domain/to the computational 
domain can have very complex shapes and can be two or three dimensions. The 
two-boundary method can map all of the boundaries of a physical spatial domain 
correctly if the two curved boundaries connect two straight lines (in the two- 
dimensional case) or flat surfaces (in the three-dimensional case). 

Application of the two-boundary method based on transfinite interpolation 
for generating grid points involves the following seven major steps (refs. 1 
and 2): 

1. Define the nature of the coordinate transformation. 

2. Select a time stretching function. 

3. Select the two boundaries of the physical spatial domain that do not 
touch each other at any point and must be mapped correctly. 

4. Describe the two boundaries selected in parametric form in coordinates 
of the transformed domain. 

5. Define curves that connect the two boundaries by using transfinite 
interpolation . 

6. Discretize the domain. 

7. Control distribution of grid points by using stretching function. 

The four-boundary method is used to algebrai cal ly generate grids for prob- 
lems in which four boundaries of the spatial domain must be mapped correctly 
from the physical spatial domain to the computational domain to the computa- 
tional domain. The four boundaries that must be mapped correctly can have com- 
plex shapes and can be two or three dimensions. The four-boundary method is an 
extension of the two-boundary method (ref. 1). The procedure for applying the 
four-boundary method is similar to the seven steps described earlier for the 
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two-boundary method. The four-boundary method involves mapping the four bound- 
aries of the spatial domain correctly to the transformed computational domain. 

In this investigation, the two-boundary method is applied to generate grid 
points for the two-dimensional step dump combustor, the two-dimensional venturi 
nozzle, the two-dimensional step dump combustor with wedge-shaped solid bodies 
embedded inside the physical spatial domain. 


ALGEBRAIC GRID GENERATION FOR A STEP DUMP COMBUSTOR 

The two-boundary method is applied to the two-dimensional axisymmetric 
step dump combustor. Figure 1 shows the physical spatial domain of the two 
dimensional axisymmetric dump combustor. Because of the symmetry of this prob- 
lem, grid points are generated for only half of the dump combustor. The physi- 
cal spatial domain consists of region A enclosed by curves 1, 2, 3 and 4 and 
region B enclosed by curves 5, 6, 7 and 8 as shown in figure 1. If the upper 
curves of the spatial domain shown in figure 1 are replaced by a single curve 
displaying a sudden step at location B, the spatial derivatives of this/curve 
forming the upper boundary of the domain will contain discontinuities at loca- 
tion B. Then the metric coefficients of the transformation needed to obtain 
solution to the conservation equation governing the problem will be discontinu- 
ous at grid points on the boundary at location B and this discontinuity on the 
boundary will propagate into the interior of the spatial domain. Therefore the 
spatial domain of interest is divided into regions A and B which are patched 
together across common boundaries as shown in figure 1. It is necessary to 
match the grids across the common boundaries and to maintain control over the 
grid spacing in these regions. 

In the following, the seven steps of the two-boundary method are applied 
to generate grid points for the step dump combustor. 


Step 1 - Define the Coordinate Transformation 

Since the two-dimensional spatial domain shown in figure 1 is nondeform- 
ing and the grid points in the x-y-t coordinate system will not move, the 
mapping of grid points in regions A and B involves the following coordinate 
transformation: 

Region A : 

1 < E _ IL 
1 < n i JL1 

(x A ,y A ,t) ^ <5,n.t> <1> 

( 2 ) 


where 


t = t<T> 
x A = x A (£,n> 


( 3 ) 



YA = yA ( S.n> 


(4) 


Region B: 


Hi < 5 < IL 
JL1 < n < JL 

<XB,yB.t)^<^,n,x) (5) 


where 


t = t ( t ) 


( 6 ) 


X B = 


(7) 


ye = ye^.n) (8) 

where x-y-t represents the coordinate system of the physical domain and 
represents a boundary-fitted coordinate system of the transformed domain 
(fig. 1). 


Step 2 - Select a Time Stretching Function 

Since there is no stretching in time for both regions A and B, the rela- 
tionship between t and x is to be taken 

t = T (9) 


Step 3 - Select Two Boundaries for Each Region 

For the spatial domain shown in figure 1, curves 1 and 2 in region A are 
chosen to correspond to coordinate lines n = 1 and n = JL1 , respectively. 
Curves 5 and 6 are chosen to correspond to n * JL1 and n = JL, respectively 
(fig. 1 ) , i .e. , 

Region A : 


*1 

= x<5,n 

= 1) = 

Xi<5> 

( 

10) 

Yl 

= y<£,n 

= 1) = 

Yl<5> 

( 

11) 

X 2 = 

; x(E,n = 

JL1 ) 

= x 2 <5> 

( 

12) 

Y z = 

y<£,n = 

JL 3 ) 

= y 2 (§) 

(' 

13) 

x 5 = 

x(E,n = 

JL1 ) 

= X 5 (E) 

(1 

14) 

y 5 = 

y<^,n = 

JL1 ) 

= T 5 <£> 

(1 

5) 
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X 6 = x<s,n = JL) = X 6 (U 
Y 6 = y(5,n = JL) = Y 6 (U 


( 16 ) 


(17) 

The two boundary method can map all of the boundaries of the physical domain 
correctly if the remaining boundaries are straight lines as in the present 
problem. As will be shown in step 4, the remaining boundaries, curves 2, 4 of 
region A and curves 7, 8 of region B, are mapped to coordinate lines 5=1, 

5 = IL and 5 = III and 5 = IL, respectively. 


Step 4 - Define the Two Boundaries Selected in Parametric Form 


Once the two boundaries have been chosen, the next problem is the repre- 
sentation of these two boundary functions in parametric form in terms of the 
parametric variable 5 . Note that since boundary lines 3 and 4 are mapped to 
coordinate line 5 = 1 and 5 = IL, respectively, the variable can only vary 
between 1 and IL ( note that x varies between A and C, see fig. 1). 

One approach is to choose the parametric variable 5 such that 


X - X »1n 
- V - X 
max min 


£, min < 5 < 5 

’mm - - 


max 


(18) 


By substituting equation (18) into equations (10) to (17), the desired paramet- 
ric equations are obtained. 


Region A : 

Xj(5) = X 2 <5> = (B - A)5 + A, 1 < 5 - < IL1 (19) 

X ] ( 5 ) = X 2 < 5 ) = (C - B)5 + B, IL 1 < 5 = *-f-| < IL (20) 

Y } ( 5 ) = D , 1 < 5 < IL (21 ) 

Y 2 (5) = E, 1 < 5 < IL (22) 

Region B : 

X 5 <£) = (C - B)5 + B, IL1 < 5 = £-f-§ < IL (23) 

Y 5 (?) = E (24) 

Xg< 5 ) = (C - B)5 + B, IL1 < 5 = < IL (25) 

Y g (£) = E + a (26) 

where A, B, C, D, E, and a are constants defined in figure 1. 
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Step 5 - Define Curves That Connect the Two Boundaries 

By using transfinite interpolation based on Hermite interpolation to 
interpolate between the boundaries selected in step 3, (i.e., curves 1, 2 and 
curves 5, 6), the connecting curves can be readily obtained. Since the bounda- 
ries curves are themselves functions and can be determined independently, equa- 
tions (3), (4) and (7), (8) can be rewritten as cubic polynomials. 

Region A : 

X A (£,n) = Xj (£)hj < n ) + X 2 (£)h 2 <n> 


3x(£ ,n 
3n 


— h 3 (n> 


+ 


a X(£,n ~ JL1 ) 
3n 


h 4<n ) 


(27) 


Y A (?,n) = Y 1 <?)h 1 (n> + Y 2 (^)h 2 ( n ) 

+ = — h 3 ( n ) + 8 y( ^'a n = JL ~ n h 4 ( n> (28) 

Region B : 

X B (5,n> = X 5 (^)h 5 ( n > + X 6 <5)h 6 ( n ) 

+ — h 7 ( n ) + 8 ~ ( ^ 1 a ' n " ~ L ~ h 8 (n) <29) 

Y B (^ >n ) = Y 5 (£)h 5 <n> + Y g Ce>h 6 (n> 

+ 9 y( ^'a n = JLn h 7 < n ) + 8 y( ^'a n = ~ h g <n> (30) 

where x a (£,n>, y a (5,n>, xg(£.n), and yB(C.n) are connecting curves: X], Yj, 

x 2 - Y 2 > x 5 > y 5 > x 6 > ancl y 6 are given by equations (19) to (26); h^(n), h 2 (n>, 
h 3 ( n ) . h 4 (n), h 5 (n), hg(n), hy(n>, and hg(n> are the Hermite polynomials 
(ref. 3). 

Applying a cubic connecting function implies that the physical grid can be 
forced to be orthogonal at the boundaries. In order for the connecting curves 
to connect curve 1 perpendicularly in the physical domain, the dot product of 

ejr (the vector tangent to curve 1 ) and In (the vector tangent to the connect- 
ing curve) at any point on curve 1 must be zero, i.e., 



or 


dX^) + + 

6l 1 + d! J 


rMLtn 

3n 


12 t + MIL lQ. 

3n 


— !] = o 


(31) 


_1 3x(^,n = 1) 1_L 3y(£ ,n = j> 

d£ 3n + d£ 3n 


( 32 ) 
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Applying this procedure will force the grid to be orthogonal at the boundary 
curve 1. Equation (32) gives the following two conditions: 


w fn n dY,<5> 

Hn - - K l (?> ~V <3 

sure i \ dX,(C) 

Y ? 8n ' - »¥*> ~k~ <3 

Control of the orthogonality for the connecting functions at other boundaries 
2 , 5 , and 6 is accomplished in a similar manner and yields the following 
conditions: 


3x(E,n = OLD 
3n 

« -K 2 <5> 

dY 2 (E> 

dE 

(35) 

3y(E,n = JL1) 

= +< 2 (C) 

dX 2 (U 

(36) 

3n 

dE 

3x(E,n = JL1) 
3n 

= -K 5 (E) 

dY 5 (0 

dE 

(37) 

3y(E,n = JL1) 
3n 

= +K 5 ( ^ 

dX 5 <C) 

d'E 

(38) 

3x(E,n = JL) 

" - K 6 ( ^ 

d V 5> 

(39) 

3n 

dE 

3y(E,n = JL) 

" +K 6 ( ^ 

dX g (E> 

(40) 

3n 

dE 


K i < f; ) , K?<C), << 5 ( 5 ), and Kg(i;) can be related to the magnitude of the boundary 
curves and are chosen by trial and error to ensure the grid lines do not over- 
lap each other at the interior of the physical domain. For the present problem 

Ki<5> = K 2 (£) = K f [Y 2 (5> - Y i ( £ ) ] (41) 

K 5 (C) = K 6 (5> = Kf[Y 6 <5> - Y 5 <C>] <42) 

Substitution of equations (33) to (42) into equations (27), (28) and (29), (30) 
gives the desired connecting curves. Given the connecting curves and paramet- 
ric boundary curves, a uniform computational grid can be mapped onto the physi- 
cal domain forming a physical grid. 


Step 6 - Discretize the Domain 

The time domain in the E-n-t coordinate system is discretized by replac- 
ing it with equally-incremented time levels, i.e., 
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n = 0, 1 , 2, . . 


x n = nAx , 


(43) 


where xi denotes time at the nth time i eve i anc j Ax j s the time-step size. 


The spatial domain in the 5-n-x coordinate system is discretized by 
replacing it with equally spaced grid points. The locations of the grid points 
are given by 

Region A : 


Region B : 


= i r 

i = 1 

, 2, . . ., 

, I LI , 


, IL 

(44) 

= j. 

3 = 1 

, 2, . 

, JL1 , 


, JL 

(45) 

Ci - i . 

i = IL1 , 

• • • i 

, IL 


(46) 

nj = j, 

j = JL1 , 


JL 


(47) 


Substitution of equations (44) and (45) into equations (27), (28) and equa- 
tions (46) to (47) into equations (29), (30) gives the grid points in the phys- 
ical spatial domain. 


Step 7 - Control Distribution of Grid Points 


The two-boundary method is applied to generate grids inside the step dump 
combustor. Boundary layer flow near the wall of the combustor is characterized 
by strong vi sci d-i nvi sci d interactions; in addition, one expects/steep velocity 
gradients at the injector location upstream of the rearward-facing step region. 
Thus, more grid points should be clustered near curves 2, 5, 6, and 8 along 
x = xi n j. Concentration of grid points near curves 2, 5, 6 is accomplished by 
replacing n in equations (27) to (30) by the following equation 


(B + 1) - (B - 1 ) C ( B + 1 ) / ( B - 1)] (1 ^ 
n n n q 

[ ( B + 1 ) / ( B - + 1 

n n 


(48) 


In the above equation, B n is a constant greater than unity. To concentrate 
more grid points near curve 8, 5 in equations (29) and (30) is replaced by the 
fol lowi ng equation 


(B + 1 ) [ ( B + 1 ) /( B - 1)] 
n n n 


( 2 ^- 1 ) 


B + 
__n 


2 { 1 + [(B + 1 )/(B - 1 )] 

n n 


(25-1 ) 


(49) 


To concentrate grid points along x = x ^ n j , replace 5 in equations (27) 
and (28) by the following expression 


si nh [ B 


«lnj 


11 


«>] 


si nh(aB^) 


(50) 
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where 


„ 1 1 n 

1 + (e * - 


a 2B^ 1n 

~ B r 

1 + (e * - 

n «ln 3 . 


In the above equations, 8^ is a constant greater than zero. By varying the 
parameter B n and B^, is different distributions of grid points can be 
obtained next to curves 2, 5, and 6 and along curve 8 and x = x} n j. 

Figure 2 shows grid points generated inside the axisymmetric step dump 
combustor with no stretching functions used and with 51 grid points spanning 
the £ direction and 21 grid points spanning the n direction. The grid 
system shown in Figure 3 was generated by using stretching functions (7) to 
(9). Approximately 50 sec of the IBM 370 CPU time is required to generate the 
grid system shown in figure 2. 


ALGEBRAIC GRID GENERATION FOR A VENTURI NOZZLE 

The two-boundary method is applied to the axisymmetric venturi nozzle. 
Because of the symmetry of this problem, grid points are generated for only 
half of the nozzle. Figure 4 shows the physical domain of the nozzle geometry. 
By following step 1 to 7 described above, a system of grid points is generated 
in the physical spatial domain. The following set of parametric boundary equa- 
tions are used (step 4) 


X] <£) = (C - A)£ + A 


(51) 


Y](£> = D 


(52) 


X 2 (£) = (C - A)£ + A 


(53) 


Y 2 C£) = 


(E 


a l } 


+ a, cos 


rr[ (C - A)£] 
8 A 


<C<^ <54. 


= (E + a 2 ) 


a 2 cos 


ir(C - A)£3 
C - B 


< £ < IL (55) 


£ = 


x - A 
c - A 


where A, B, C, D, E, a ] , and a 2 are constants defined in figure 4. 

Figure 5 shows grid points generated inside the axisymmetric nozzle with 
no stretching function used and with 51 points spanning the £ direction and 
21 grid points spanning the n direction. Approximately 50 sec of the IBM 370 
CPU time is required to generate the grid system shown in figure 5. The non- 
uniformly distributed grid system shown in figure 6 was generated with the use 
of stretching function. Comparing figure 5 and 6 indicates that the stretching 
functions employed can cluster grid points near the wall of the nozzle. 


x - A 

c - A 
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ALGEBRAIC GRID GENERATION FOR A STEP DUMP COMBUSTOR 


WITH EMBEDDED SOLID BODIES 

To evaluate the feasibility of the two-boundary method, a two-dimensional 
axisymmetric step dump combustor with geometrical ly-complex embedded solid bod- 
ies is considered. The flow field domain is divided into five subdomains as 
shown in figure 7, each covering only part of the whole field. Grid points are 
generated inside each of the subdomain using the two-boundary method. The grid 
systems of these subdomains are interfaced with each other in such a way that 
the grid lines are continuous at the common boundaries. 

Figure 7 shows the physical domain of the step dump combustor geometry 
together with the partitioned subdomains. The two-boundary method is applied 
to each of regions A, B, C, D, and E. The boundary curves are defined as 
described below. 

Region A : 

IL4 < ? < IL 
1 < n < JL1 

Curve 1 : 


X } <V = E + (F - E)(jrf-f) 

Y, (?) = P 

Curve 2: 

X 2 (?> = E + (F - E)(jU-|) 

y 2 <?) = 0 

Region B : 

1 < ? < IL 
JL1 < r, < JL2 

Curve 5: 

X 5 (?) = A + (B - A)(|-5-J), 1 < ? < IL1 

= B + (C - IL1 < ? < IL2 

= C + (D - C > Z“§) . IL2 < ? < IL3 


(55) 

(56) 

(57) 

(58) 


(59) 

(60) 
(61) 


10 



= D + (E - D>(|-f-§), IL3 < 5 < IL4 (62) 

- E + <F - E>(|-5-|), IL4 < ? < IL (63) 

Y 5 (E> = P, 1 < E, < IL2 (64) 

- P + (Q - P>(g-fr). IL2 < E < IL3 (65) 

= Q, IL3 < E < IL (66) 

Curve 6: 

Xg(£) = X 5 <C), 1 < K < IL1 (67) 

= X 5 <?), IL1 < Z <_ IL2 (68) 

= X 5 (E>, IL2 < i < IL3 (69) 

= Xg<E), IL3 < E < IL4 (70) 

= X 5 <E), IL4 < £ < IL (71) 

Yg(E> = S, 1 < E < IL2 (72) 

= S - (S - IL2 < E < IL3 (73) 

= R, IL3 < E < IL (74) 


Region C : 

IL4 < ^ < IL 
JL2 < n < JL3 

Curve 9: 

x 9 (5) = E + (F - E)(f^-|) 

Y 9 (E) - E + (F - E)(^-|) 

X 10 (E) = X g (?) 

V 5) = T 
1 1 


( 78 ) 



Region D : 


Curve 13: 


Curve 14: 


Region E : 


Curve 17: 


1 < £ < IL 
JL3 < n < JL4 


X 13 ( C) = X 5 <C), 
- x 5 <e> * 


= x 5 <o 


= x 5 <5). 


- 


1 < 5 < IL1 

IL1 < 5 < IL2 

IL2 < K < IL3 

IL3 < K < IL4 

IL4 < i < IL 


Y 13 <U - S, 1 < C < IL2 


= S + (T - S>(j^-§), IL2 < i < IL3 
= T, IL3 < 5 < IL 


X 14 <5) = x 5 <?) 

= X 5 (P 


= v 5> * 


- V 5> 


- 


y 14 (C) = u. 


1 < C < IL1 

IL1 < i < IL2 

IL2 < C < IL3 

IL3 < £, < IL4 

IL4 < £ < IL 

1 < ^ < IL2 


IL1 < £<IL 
JL4 n <JL 


(79) 

(80) 
(81) 
(82) 

(83) 

(84) 

(85) 

( 86 ) 

(87) 

( 88 ) 

(89) 

(90) 

(91 ) 

(92) 


X 17 (5) = X 5 (?), 1 < ? < IL2 


(93) 
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(94) 


= X 5 U>, IL2 < i < IL3 


= Xg(5), IL3 < £ < IL4 (95) 

= X 5 (5), IL4 < (■ < IL (96) 

Y iy (5> = U, IL1 < C < IL (97) 

Curve 18: 

X 18 (£> = X s ( 5 ). 1 < C < IL2 (98) 

= X 5 (C), IL2 < i < IL3 (99) 

= X 5 (£), IL3 < i < IL4 (100) 

= , IL4 < ^ < IL (101 ) 

Y la <C> = V, IL1 < 5 < IL (102) 


where the constants appearing in the above equations are defined in figure 7. 
The functional form for the connecting curves can be readily obtained by using 
Hermite interpolation to interpolate between the two boundaries selected for 
each of the subdomain. 


Region A : 

X A (£,n> = X 1 (C)h 1 ( n ) + X 2 (£)h 2 (n> 


+ 8 -- ^ - ]) h 3 ( n ) + ax( ^a n ~" L1> h 4 (r i> (103) 
Y A (C,n) = Y ] (5)h ] (p) + Y 2 (C)h 2 ( n ) 

+ 8 — ^ n - - - 1 - h-( n > + 8 - ( ^ = — L1) h,<n> (104) 

3n 3 on 4 

Region B : 

X B (£,n> = X 5 (^)h 5 ( n ) + X 6 (?)h 6 (n> 

+ 9x< 5 « H— - JLl) h ( ) + 9x(^,n = J L2) h ( } (105) 

on 7 on 8 

Y B <5, n ) = Y 5 <?)h 5 (n) + X g (e)h 6 (n) 

, • JLi) ( , , 3x( . ;,n - J L2) . <l06) 

dp 7 dp o 
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Region C : 


X c <5.n) = X 9 (C)h g ( n ) + X 10 (C)h 10 ( n ) 


. 3x(C,n = JL2) , ^ 3x(^,n = 313) , , x 

aTT h n ( n> + h 10 (n> 


3n 12 


Y c U,n) = Y 9 <C>h 9 <n) + Y ]() (C)h 10 <n) 


3 x(C, n = JL 2 ) . , x 3 x(£ ,n = JL 3 ) _ , x 

* h u <n) + — aq h i 2 < ’ 1> 


Region D : 

x D <C,n> = X 13 (§)h 13 ( n ) + X 14 <c>h 14 ( n ) 


3L4) 


3n 15 

Y D (5,n> = Y 13 (§)h 13 (n) + Y 14 (C)h ]4 ( n ) 


«■ JL3> hic(n) + - 8 -- ^ n — u;T ' h 16 ( n ) 


Region E : 

x E (^,n) = X 17 (£)h 17 ( n > + X 18 (5)h 18 ( n ) 

, 3x(£,n = JL4) . . x 3x(£, n = JL) L , v 

+ ~3^ h 19 (n) + ^3n h 20 n 

Y E <?,n) = Y 17 (e>h 17 ( n ) + Y 18 <c)h 18 ( n ) 


( 107 ) 


(108) 


(109) 


( 110 ) 


(111) 


+ 


3x(£,o = JL4) 
3n 


hi g(p) 


3x(g, n = JL) 
+ 3n 


h2g(nl 


(112) 


The Hermite polynomials are described previously derivative terms appear- 
ing in equations (103) to (112) are chosen so that the connecting curves per- 
pendicularly intersect the boundaries selected above. Grid points inside this 
step dump combustor configuration are generated by using equations (103) 
to (112). Figure 8 shows grid points generated inside the spatial domain of 
this step dump cumbustor with 51 grid points spanning the £ direction and 
21 grid points spanning the n direction. Approximately 50 sec of the IBM 
CPU time is required to generate the grid system shown in figure 8. Figure 9 
shows the grid system generated with the use of stretching functions. 
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ALGEBRAIC GRID GENERATION FOR A RECTANGULAR DUCT WITH ARBITRARY 


CROSS SECTIONS AND EMBEDDED SOLID BODIES 

In this section, the four-boundary method (ref. 1) is applied to generate 
grid points for an axi symmetr i c rectangular duct with arbitrary cross sections 
and embedded solid bodies as shown in figure 10. 

The grid is case in three dimensions by treating the three-dimensional 
spatial domain as a stack of two-dimensional spatial domains. The three- 
dimensional grid is generated by generating two-dimensional grids at equally 
incremented locations along the z-axis in the x-y-z-t coordinate system. 
Because of the symmetry of this problem, grid points are generated for only 
half of the duct shown in figure 10. The three-dimensional grid is obtained 
by rotating the upper half description about the axis of symmetry. Figure 11 
shows the physical domain of the two-dimensional spatial domain. The spatial 
domain of interest is divided into regions A and B. The four-boundary method 
is applied to each of regions A and B by considering curves 1, 2, 3, and 4 as 
boundary curves of region A and curves 5, 6, 7, and 8 as boundary curves of 
region B. 

The boundaries are defined as described below. 

Region A : 


1 £ S <. IL 
1 < n < JL1 


Curve 1 : 


x, (£,n = 1) = x 1 <c> = [fiTTT <R + + R_>] - R 
y j (£ » n = 1 ) = Y 1 (£) = R tan 


Curve 2: 


x 2 (E,n = JL1) = X 2 (£) 


\i -_.L 

LI L - 1 


(R + R ) 
+ - . 


R 


(113) 

(114) 


(115) 


y 2 (£,n = J LI ) = Y 2 <t;> = R tan a 2 


(116) 


Curve 3: 


y 2 (£ = I,n) = Y 3 (n> 
Curve 4: 


x 3 (£ = 1 ,n) = X 3 (n> = -R_ 


(R tan a 2 - R tan ) 


+ R tan a, 


(117) 

(118) 


x 4 (£ = I L , n ) = X 4 ( n ) = R + 


(119) 
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( 120 ) 


y 4 (£ - IL, n > = Y 4 < n > = [ juT T (R tan a 2 ~ R tdn a 1 } ] + R tan a l 
Region B : 


1 < 5 < IL 
JL1 < n < JL 


Curve 5: 


x 5 (c,n = jli> = x 5 <e> = [frrr (R + + R _ } 


R 


= JL1 ) = Y^(£) = R tan 


Curve 6: 


x 6 (C,n = JL) = X 6 (V = [f^TT (R + + R ->] ' R - 
y 6 (^,n = JL) = Yg(5> = R tan <x 4 


Curve 7: 


x y (5 = 1 ,n) « X 7 <n> = -R 


y 7 ( 4 

= l,n> 

II 

-< 

/^v 

J3 

ii 

fn - 
Ljl - 

JL1 
• JL1 

(R 

tan 

a 4 

R tan a^) 

+ R 

tan 

a 3 

Curve 

8 : 











OO 

sTr\ 

II 

IL,n) 

ll 

X 

00 

.3 

II 

rn_=_ 

Ljl - 

JL1 
• JL1 

(R 

tan 

a 4 " 

R tan a.j) 

+ R 

tan 

a 3 




y 8 <{ 

= IL 

,n) 

= Y 

8 ( n) 

= R 





( 121 ) 

( 122 ) 

(123) 

(124) 

(125) 

(126) 

(127) 

(128) 


where R, R +t R_ a ] , « 2 , 013 , and a .4 are defined in figure 11. 

By using Hermite interpolation, cubic polynominals are used to generate 
connecting curves for the boundary curves. The connecting functions are 


Region A : 


X A (£,n> = X (12) <4,n) + X VJtw (4,n) - x"^'(4 = 1.n>h 5 (4> - X VI ^'(4 = IL,n>h 6 <4> 


,(34) 


,( 12 ) 


( 12 ), 


a ( 12 ), _ , . 
3x ( £ s 1 > q ) 

94 


„ (12) lt . Tl s 

h 7 (?) - — hg(U 


94 


(129) 
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Region B : 


Y A <5, n > = Y (12> (C.n> + Y (34) (^,n) - y ( 12) (£ = l,n>h s <5> - Y (12) (^ * IL,n>h 6 <5> 

_ 8y< 12> g| = • ' n - h 7 <C) - ax<12)( ||~ IL > - ri - hg(5) (130) 

where 

x (12) (s, n > = X 1 <?)h 1 <n) + X 2 <5)h 2 <n> 

+ ax( ^ = h 3 < n ) + 9x( ^^° - 3L --- h 4 <n> (131) 

X ( 34 ) (£,n) = X 3 ( n )h 5 (^) + X 4 ( n )h g (C) 

+ 3x( ^= 1in - - h ? <£) + 9x( ^ ^ -^ n ~ h g <5) (132) 

Y ( 12 ) (S >n > = Y 1 ( f; ) h 1 (g) + Y 2 (^)h 2 ( n ) 

+ — h 3 (n> + 8y< ^ , a n °"""~ h 4 < n) (133) 

Y (34) (^, n ) = Y 3 ( n >h 5 (ip + Y 4 < n )h 6 <C> 

+ hy(^) + 8y( ^ hg<S> (134) 

h]<5), h 2 (i), h 3 <?), h 4 <C(, h 5 <£), h 6 (£), h 7 (U, and h 8 (£) are Hermite polyno- 
mials and are described previously. 

The mapping between the physical domain and the transformed computational 
domain is obtained by using equations (113) to (134). Figure 12 shows the 
transformed domain in three dimension. Figure 10 shows grid points generated 
inside the axisymmetric rectangular duct with arbitrary cross section and solid 
bodies embedded in the interior physical domain with 41, 21, 31 grid points 
spanning the £, n. and (, directions, respectively. Approximately 100 sec of 
the IBM 370 CPU time is required to generate the gird system shown in figure 10. 


CONCLUSION 

Algebraic grid generation method called the two-boundary method has been 
used to generate grid points inside the spatial domains of an axisymmetric step 
for dump combustor, an axisymmetric venturi nozzle and a step dump combustor 
with wedge-shaped solid bodies embedded inside the physical spatial domain. 
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Grids are generated for the three-dimens ional rectangular duct with arbitrary 
cross sections and with solid bodies embedded inside the spatial domain using 
the four-boundary method. Grid generation for these complex arbitrary geome- 
tries can be obtained at a reasonably fast computer rate. The successful 
generation of grids for the cases studied demonstrates the usefulness and via- 
bility of the two-boundary and four-boundary algebraic grid generation methods 


REFERENCES 

1. Shih, T.I-P.: Fi ni te-Defference Methods in Computational Fluid Dynamics. 

To be published by Prentice Hall, Englewood Cliffs, NJ. 

2. Thompson, J.F.; Warsi, Z.U.A.; and Mastin, C.W.: Boundary Fitted 

Coordinate Systems for Numerical Solution of Partial Differential 
Equations - A Review. J. Comput. Phys., vol . 37, July 1982, pp. 1-108. 

3. Yang, S.-L.; and Shih, T.I-P.: An Algebraic Grid Generation Technique for 

Time-Varying Two-Dimensional Spatial Domains. Int. J. Numer. Meth. 

Fluids, vol. 6, May, 1986, pp. 291-304. 

4. Thompson, J.F., ed.: Numerical Grid Generation. North Holland, 1982. 

5. Smith, R.E.: Two-Boundary Grid Generation for the Solution of the 

Three-Dimensional Compressible Navier-Stokes Equations. NASA TM-83123, 

1981 . 

6. Vinokur, M.; and Lombard, C.K.: Algebraic Grid Generation with Corner 

Singularities. Advances in Grid Generation, vol. 5, K.N. Ghia and Vrmila 
Ghia, eds., ASME, 1983, pp. 99-106. 

7. Dolce, J.S.; Shih, T.I-P.; and Roan, V.P.: One-Dimensional Stretching 

Functions from a Variational Principle. CFDL Report No. 8801, Department 
of Mechanical Engineering, University of Florida, 1988. 

8. Vinokur, M.: On One-Dimensional Stretching Functions for Finite-Difference 

Calculations. J. Comput. Phys., vol. 50, May, 1983, pp. 215-234. 

9. Anderson, D.A.; Tannehill, J.C.; and Pletcher, R.H.: Computational Fluid 

Mechanics and Heat Transfer. Hemisphere Publishing Corp., Washington, 

1984. 

10. Mastin, C.W.; and Thompson, J.F.: Elliptic Systems and Numerical 

Transformations. J. Math. Anal. Appl . , vol. 62, 1978, pp. 52-62. 

11. Thompson, J.F.; Thames, F.C.; and Mastin, C.W.: Boundary-Fitted 

Curvilinear Coordinate Systems for Solution or Partial Differential 
Equations on Field Containing Any Number of Arbitrary Two-Dimensional 
Bodies. NASA CR-2729, 1977. 


18 



12. Thames, F.C., et al . : Numerical Solutions for Viscous and Potential Flow 

About Arbitrary Two-Dimensional Bodies Using Body-Fitted Coordinate 
Systems. J. Comput. Phys., vol . 24, July 1977, pp. 245-273. 

13. Hirsh, R.S.: Higher Order Accurate Difference Solutions of Fluid 

Mechanics Problems by a Compact Differencing Technique. J. Comput. Phys., 
vol. 19, Sept. 1975, pp. 90-109. 

14. Smith, G.D.: Numerical Solution of Partial Difference Equations: Finite 

Difference Methods. Second Ed., Clarendon Press, Oxford, 1978. 

15. Steger, J.L.: Implicit Finite-Difference Simulation of Flow About 

Arbitrary Two-Dimensional Geometries. AIAA J., vol. 16, no. 7, July, 

1978, pp. 679-686. 




FIGURE 1. - (A) THE PHYSICAL SPATIAL DOMAIN. CURVE 1 IS DESCRIBED BY = D WHEN X £ 
[A,C3. CURVE 2 IS DESCRIBED BY Y 2 = E WHEN X £ [A,CL CURVE 5 IS DESCRIBED BY Y 5 - E 
WHEN X £ IB, CL CURVE 6 IS DESCRIBED BY Y & = E + a WHEN x £ [B,CL A, B, C, D, E, 

AND a ARE CONSTANTS. (B) TRANSFORMED COMPUTATIONAL DOMAIN. 
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FIGURE 4. - (A) PHYSICAL SPATIAL DOMAIN. CURVE 1 IS DESCRIBED BY Y, = D WHEN x £ IA,CI. 
CURVE 2 IS DESCRIBED BY Y 2 = (E + a,) + a, cos [ *^ X _~^ A) ] WHEN x E IA,B) AND BY Y 2 = 
(E + a 2 ) - a 2 cos [^ x .' b B) ] wh en x £ tB.ci. a. b. c, d. e. a v and a 2 are constants. 
(B) TRANSFORMED COMPUTATIONAL DOMAIN. 



FIGURE 5. - GRID FOR VENTURIE NOZZLE. 


FIGURE 6. - GRID FOR VENTURIE NOZZLE, STRETCHING FUNCTION USED. 



FIGURE 7. - (A) PHYSICAL SPATIAL DOMAIN. (B) TRANSFORMED COMPUTATIONAL DOMAIN. 


FIGURE 8. - GRID FOR STEP DUMP COMBUSTOR WITH EMBEDDED SOLID 
BODIES. 




FIGURE 9. - GRID FOR STEP DUMP COMBUSTOR WITH EMBEDDED SOLID BODIES, 
STRETCHING FUNCTION USED. 



FIGURE 10. - GRID FOR RECTANGULAR DUCT GEOMETRY WITH EMBEDDED SOLID 


(B) 


FIGURE 11. 
DOMAIN. 


"Wl 

(A) 

- (A) TWO-DIMENSIONAL PHYSICAL SPATIAL DOMAIN. (B) TWO-DIMENSIONAL COMPUTATIONAL 



FIGURE 12. - THREE-DIMENSIONAL TRANSFORMED COMPUTATIONAL 
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